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SUMMARY 


A tapered composite laminate subjected to tension load was analyzed 
using the finite element method. The { [0^/(±45) ]/^ [ (±45) ^]/[0/(±45)/0] 

glass/epoxy laminate has a (±45) ^ group of plies dropped in three distinct 

steps, each 20 ply- thichnesses apart, thus forming a taper angle of 5.71 
degrees. Steep gradients of interlaminar normal and shear stress on a 
potential delamination interface suggest the existence of stress 
singularities at the points of material and geometric discontinuities 
created by the internal plydrops. The delamination was assumed to initiate 
at the thin end of the taper on the -45/+45 interface indicated by the 

arrow in the laminate layup and the delamination growth was simulated in 

both directions, i.e. along the taper and into the thin region. The total 
strain-energy-release rate, G, and the mode I and mode II components of G, 
were computed at the delamination tips using the virtual crack closure 
technique. In addition, G was calculated from a global energy balance 
method. The strain-energy-release rate for a delamination growing into the 
thin laminate consisted predominantly of mode I (opening) component. For a 
delamination growing along the tapered region, the strain-energy-release 
rate was initially all mode I, but the proportion of mode I decreased with 

increase in delamination size until eventually total G was all mode II. 

The total G for both delamination tips increased with increase in 
delamination size, indicating that a delamination initiating at the end of 
the taper will grow unstably along the taper and into the thin laminate 
simultaneously. 


INTRODUCTION 


Composite rotor hubs are currently being designed and manufactured that 
are hingeless and bearingless to reduce weight, drag, and the number of 
parts in the hub. Such a design would involve tapering the laminate by 
dropping some plies in the flexure region of the hub. The plydrop in the 
laminate creates geometric and material discontinuities that create large 
interlaminar stresses and initiate delaminations. Therefore, there is a need 
to analyze tapered laminates with ply drops to understand their failure 
mechanisms . However, only a limited amount of literature is available on 
tapered laminates . 

Adams et. al. [1J analyzed a [0^/(±45)^/90^] graphite/epoxy laminate 

in which two zero degree plies were dropped. The effect of compressive load, 
moisture, and temperature due to the presence of the plydrop was studied 
using a 3-D finite element analyses with nonlinear orthotropic response. 
They concluded that all the interlaminar stresses induced by a 0 degree ply 
drop-off anywhere in the laminate were negligible compared to the in-plane 
stresses. However, they did not account for the low interlaminar strength of 
the composite compared to the in-plane strength. 

Cannon [2] conducted experiments on graphite/epoxy tapered laminates 
from the [±45/0]^ and [±15/0]^ families, subjected to tension load. For most 

laminates the failure mode and the failure stress were similar to that of 
the untapered specimen at the thin (dropped) end of the laminate. An 
analysis based on the minimization of total potential energy which accounted 
for the effect of eccentricity due to the plydrop was used to predict the 
in-plane failure stresses in unsymmetric laminates. The tests on 
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[±45/0/(±45/0)^] , where D denotes dropped plies, showed that dropping a 

number of pi ies lumped together can change the initial damage from in-plane 
failure to delamination. 

Kemp and Johnson [3] analyzed a tapered beam having a single plydrop 
using the finite element method. Symmetric and unsymmetric laminates were 
modeled as a generalized plane deformation problem subjected to a uniform 
strain in the longitudinal direction. The layups considered were 
(±45/0/90/0^/90/0/145)^ and (0/90/145/0^/145/90/0)^ where n, the number of 

dropped zero degree plies, was chosen to be 1, 2, or 3. Failure strains were 
calculated corresponding to resin failure at the dropped plies, based on a 
maximum principal stress criterion. Alternatively the failure strains were 
obtained for intralamina failure in tension and compression, using the Tsai- 
Wu criterion. The first failure event in tension or compression was 
predicted to occur in the resin. 

Although the stress dis tr ibutions in the laminate help to identify the 
highly stressed critical areas, maximum stress or strain criteria cannot be 
used to predict delamination onset and growth if the stresses are singular. 
However, interlaminar fracture toughness, which is generic to a given 
composite material, can be used to predict the loads corresponding to the 
onset and propagation of delamination [4,5,6]. For example, delamination 
growth can be predicted from the mode I and mode II components of the 
strain-energy-release rate under static loading and from the total strain- 
energy- release rate for fatigue loading [5,7]. Therefore, the purpose of 
this paper is to study the interlaminar stress distributions in a tapered 
beam subjected to tension loads and to determine the strain-energy-release 
rate for delamination growth that may occur due to the presence of plydrops. 
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A typical stacking sequence used in a helicopter hub Is 
l [0g]/[ (±45) j]/[ (± 45 ) 2 ! ) g . The laminate considered here is a { [ 0_,/(±45) ] / 

[ (±45) ^ ] / [ 0/(±45 ) /0 ] } ^ tapered laminate. This laminate has the same number 

of 0° and ±45° plies as the hub but a somewhat different stacking sequence. 
The (±45)^ plies are dropped in three steps, 20 ply thicknesses apart. The 

dropped plies result in a taper angle of 5.71°. The laminate was analyzed 
using a two-dimensional finite -element analysis. The interlaminar normal and 
shear stress distributions along the taper interface, indicated by an arrow 
in the above layup, are presented. Delaminations are assumed to initiate at 
the point of highest interlaminar stress along this interface. The mode I, 
mode I_I and total strain-energy-release rates for various delamination 
lengths are presented. These results were used to hypothesize the stability 
of delamination growth under static and fatigue loading. 


NOMENCLATURE 


a 

b 

E ll, E 22, E 33 
G 


G I' G II, G III 


delamination length along taper 
delamination length in the thin region 
Young's moduli 

total strain-energy-release rate 

mode I, mode II, and mode III components of 


G 


12 * 


G 13 ’ G 23 
h 


st rain -energy -re lease rate , respectively 
shear moduli 

ply thickness 

total load per unit width on symmetric half laminate 
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distance along delaminations 
Cartesian coordinates 
uniform tension load per unit area 

interlaminar normal stress 

interlaminar shear stress 

Poisson's ratios 

ANALYSIS 

Specimen Configuration and Loading 

Figure 1 shows the tapered laminate that was analyzed. The stacking 
sequence Is { [0 ? /(±45) ]/ [ (±45) 3 ] / [ 0/(±45)/0 ] } . The (0 ? /±45) ply group in 

the laminate of Fig.l forms the belt area, and the (0/145/0) laminate In 

the center forms the core. The transition from the thick region at the left 
to the thin region at the right is achieved by dropping the group of (±45) ^ 

plies in three distinct steps, each 20-ply thicknesses apart. The shaded 

regions shown in Fig.l are the resin pockets formed at the ends of the ±45 
degree plies that are terminated. In similar laminates, delaminations have 
been observed at the interface indicated by the arrow in the layup above. 
Therefore, the delaminations are assumed to grow along the interface ABCD in 
Figure la. A typical delamination is shown in Fig. lb. The delamination Is 

assumed to form at point C, and grows into the tapered region (tip I) and 

into the thin region (tip H) . 

The tapered laminate was assumed to be made of S2/SP250 glass/epoxy and 
to be subjected to a uniform load at the thick end (X«0) . Examination of the 


x 

X, Y, Z 


n 


nt 
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results indicates that the displacements are uniform in the neighborhood of 
X~60h, Thus, the uniform load condition at X=0 is equivalent to a uniform 
displacement condition. A fixed grip condition was assumed at the thin end. 
The material properties used in the analysis are given in Table 1. The in- 
plane properties for a unidirectional ply (e.g; E^, ^22* ^12’ ^12 ^ are 

similar to those used in reference 7. The out-of -plane properties (G^, ^13 ' 
^23* ^23 ^ were assumed to be identical to the in-plane properties, and 

was assumed equal to ^2* 

Finite Element Model 

A 3-D finite element analysis of the laminate is desirable, but such 
analyses are complex. Simple 2-D models, which do not account for the free 
edges, usually provide insight that can be used in 3-D analyses. Thus, as a 
first step, 2-D plane-strain analyses were performed in this study. 

Furthermore, the stacking sequence considered here contains only 0 degree 
and ±45 degree plies. With the absence of the 90 degree plies, the 
interlaminar Poisson mismatch between plies that causes edge delaminations 
was not considered to be significant [5]. Therefore, a two dimensional 
finite -element analysis is expected to be reasonably accurate for this 
laminate . 

A two dimensional finite element model was developed utilizing the 
symmetry of the laminate about the X-axis. The model had 7610 nodes and 2382 
eight-noded, isoparametric, parabolic elements as shown in Fig. 2a. A 

refined mesh was used near plydrop points (B, E and F and C in Figure lb) to 

capture the local influence of these geometric discontinuities and the 
corresponding s tresses . The smallest element size used in the model was 
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equal to one-quarter of the ply thickness. These small elements were 
provided near the plydrops on line BC, at the transition point B from the 
thick region to the tapered region, and the transition point C from the 
tapered region to the thin region. The element size immediately below line 
BC varied in the Z-direction due to the change in the resin thickness from 
two to zero ply thicknesses in the three resin pockets. Collapsed eight- 
noded elements were used at locations E, F, and C in the resin pockets. 
Figure 2b shows local mesh detail at location E. A similar pattern was used 
at points F and C. 

The nodes at the end of the thin region (at X - 180h in Figure la) of 
the laminate were constrained in both X- and Z- directions. A uniform 
tension per unit area, cr Q , (assuming a unit width in the Y-direction) was 

applied along the X«0 line of the model. Plane strain conditions were used 
in the analysis. 

To facilitate modeling delaminations along ABCD, duplicate nodes were 
created in the model all along lines AB , BC, and CD. Multi-point constraints 
were imposed for the duplicate nodes. Different size delaminations were 
simulated by relaxing the multi -point constraints for the appropriate nodes 
along lines BC and CD. Note that two delaminations are assumed to maintain 
symmetry about the X-axis. 

The material directions of plies in the laminate are oriented at an 
angle relative to the global coordinate system of the analysis. The material 
stress - strain relations for these plies were transformed to obtain the 
stress - strain relations in the global system. Appendix A presents the 
details of the transformations used. 
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Confutation of Strain -Energy -Release Rate 

The virtual crack closure technique (VCCT) was used to obtain the 
strain-energy- release rate components, mode I, and mode II, based on the 
local forces at and ahead of, and the relative displacements behind, the 
delamination tip. These two components were calculated using the following 
equations (see Fig. 3). 


2 A 


t F ni ( 


- v> 


+ F 


n J 


( v, 


m 


v> J 


(la) 


G II - ' 2"! [ F ti ( U k ' V> + F tj < U m ' V> ] (lb) 

where A Is the element size, F ni and F ti are the normal (n) and tangential 

(t) , forces, respectively, at node i, and (v^ - v^/) and (u^ - u^/) are 

the relative opening and sliding displacements, respectively, at node k 
(see Fig. 3). Forces at node j and relative displacements at nodes m and 

m' are defined similarly. Equations 1 are similar to those given in 
references 8 and 9. The total strain-energy-release rate, G, at the 
delamination tip was calculated as 

G = Gj + G n (2) 

The mode III component of G was identically zero because plane strain 
conditions were assumed in the analyses. 
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Alternatively, the global energy change of the laminate due to 
delamination growth can also be used to calculate the total strain- energy- 
release rate, G. The strain energy of the laminate, U, can be conveniently 
computed as U- 1/2 (S u^) where f^ and u^ are the nodal forces and 

corresponding nodal displacements, respectively, for all nodes i on the line 
X=0 in Figure la. The strain-energy-release rate for successive delamination 
growth was calculated as 


dW _dU 

dA * dA 


(3a) 


where dW/dA and dU/dA are the rate of change of work and strain energy, 
respectively, with change in delamination area. In the finite-element 
analysis, Equation (3a) can be computed as 

G - < Vdx - V/ dx (3b) 


where U , and U are the strain energies for delamination lengths x+dx 
x+ax x 

and dx, respectively. The value of G thus calculated is considered to be the 
strain-energy-release rate at (x+ dx/2) , which is located at the center of 
the interval . 


RESULTS AND DISCUSSION 

First, the interlaminar stress distributions along the interface ABCD 
are presented. Next, the strain- energy- release rate variations for various 
size delaminations assumed along the interface line ABCD are shown. Finally, 
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the peak values of the total strain-energy-release rate and the mode I 
component values are presented and their significance discussed. 

Interlaminar Stresses 

Figure 4 shows the normalized interlaminar normal stress, (a^/a Q ) , 

along lines AB, BC, and CD in the laminate. Stresses were calculated in the 

local coordinate system, normal to the interface ABCD. The interlaminar 

normal stress shows peaks near the points of geometric and material 

discontinuity i.e. at points B, E, F, and C. The largest tensile value of 

the o distribution occurred at the transition point C. At the plydrops, 
n 

points B, E, F, the stresses changed from a high compressive value 
immediately to the left of the plydrop to a high tensile value Immediately 
to the right of plydrop. The variation of normalized interlaminar shear 
stress, ( r n t/ a 0 ^ 9 a ^ on S the same interfaces AB, BC, and CD Is shown in 

Figure 5. The shear stress also shows peaks at points B, E, F and C. 

These sudden changes in the normal and shear stress distributions at 
points B, E, F, and C are not unexpected. At these points, the material 
stiffness is different in different directions (see Fig. 6). Therefore, at 
points B, E, F, and C, stress singularities probably exist [10]. 

In order to investigate if this is true, a two-dimensional finite- 
element analysis of a homogeneous tapered laminate was performed with the 
same model as in Fig. 2. The tapered laminate was assumed to be of an 

isotropic material. The normalized interlaminar normal stress (a^/a^) 

distribution along the line ABCD is presented in Fig. 7. At points E and F, 
the stiffness is same in different directions irrespective of how these 
points are approached. Thus, no sudden changes in stress distribution exist 
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at these points. The normals to lines AB and BC at point B have different 
directions. Similarly, normals to lines BC and CD are different at point C. 
Thus, except for very small discontinuities at these points, the stress 
distribution all along ABCD is smooth. This confirms that the sharp changes 
in stresses observed in Figures 4 and 5 are solely due to material 
discontinuities at the points B, E, F, and C. 

Strain-Energy-Release Rate Analysis 

Delamination growth in a laminated composite structure may be predicted 
from the mode I, and mode II components of the strain-energy release rate 
under static loading and from the total strain- energy- release rate for 
fatigue loading [5,7], The computation and the use of the strain-energy- 
release rate in delamination prediction for the tapered laminate are 
discussed below. 

As seen in Figure 4, point C has the highest value of interlaminar 
normal stress, a^, compared to any other location on the interface line 

ABCD. Therefore, the delamination was assumed to initiate at this point. 
Delamination lengths a and b (see Fig. lb) were assumed within the 
tapered region along CS and in the thin laminate along CD, respectively. The 
strain- energy- release rate values G, Gj, and G^ were computed at each 

delamination tip using the finite element analysis and Equations 1-3 for 
various values of a and b . 

The total strain-energy-release rates were calculated using two 
different methods; VCCT (equation 2) and from global energy change (equation 

3). These G values normalized by N*/h, (where N y is defined as the product 
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of uniform tension stress and half the laminate depth at X * 0, and h 

is the ply thickness), are plotted for comparison in Figure 8. For this 
case, no delamination was assumed along the taper, CB, (i.e. a-0) and the 
values of G were obtained for various values of delamination lengths, b, 
along CD in the thin region of the laminate. Excellent agreement between the 
q values computed by the two methods was obtained. Similar agreement was 
found for all the cases studied. The G values obtained by using equation 3 
are presented in this paper because more data points were available for this 
computation and values of the individual modes, Gj and G^j, were taken from 

the VCCT calculation. 

Figure 9 presents a composite of G distributions for delamination 
growth in the thin and thick regions. In the right hand portion of the 
figure, the G values for the delamination tip at "H" were plotted against 
b/h for a fixed value of a/h. Similarly, the left hand portion of Figure 9 
presents the G values for the delamination tip "I" plotted against a/h for a 
fixed value of b/h. 

Referring to the delamination tip at "H" on the right side of Figure 9 
(where a/h is held constant and b/h varies) , the G initially increases 
rapidly with b/h as the delamination grows into the thin laminate along line 
CD. For a/h- 0, 6, and 12, the G attains a peak value and drops slightly 
with further delamination growth. This drop decreases with increasing a/h 
and does not occur for a/h- 20 and 24 in the range of b/h considered. These 
values of G are given in Table 2. 

In a complementary situation shown on the left side of Figure 9, the G 
values at the delamination tip I were plotted against a/h for various values 
of b/h. The total strain-energy- release rate increases initially, and then 
is relatively constant, or drops slightly, before approaching the plydrop. 
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In the proximity of the dropped plies, however, G values increase rapidly 
and attain peak values at the plydrops (a-=20h and 40h) . These G values are 
given in Table 3 . 

The results of Figure 9 suggest that a delamination initiating at point 
C will grow in an unstable manner simultaneously along the tapered Interface 
CB as well as in the thin laminate along CD. This can be explained as 
follows. Consider a small delamination initiating at point C in Fig. 9. If G 
for the left and right delamination tips control growth along the thin and 
tapered regions, respectively, then delaminations would arrest after they 
had grown to the peak values. However, growth of one delamination tip will 
increase the G for the other tip, causing growth in that direction also. 
Hence, as soon as a stable situation occurs on one delamination tip, it will 
increase the G for the other tip causing further growth. Hence, as the peaks 
in the G values on the left and right sides of Figure 9 increase 
monotonically with increasing a/h and b/h, a delamination initiating at 
point C will grow unstably in both directions simultaneously. 

As discussed above, the values of G shown in Fig. 9 reached peak values 
for delamination growth In either direction. For delamination growth along 
the thin side CD, the plots of G for delamination tip "H" vs. b/h are 
similar to those obtained for edge delamination, where G is initially zero 
at b/h=0 and increases to a plateau at some distance, usually b/h- 2 to 3 
[5-7]. The distance at which G reaches a plateau for the edge case may vary 
with the interface analyzed [5]. Similarly, the distance at which G at tip 
"H" vs. b/h reaches a peak varies with a/h. The value of b/h to reach this 
peak, however, is of little consequence since it is assumed that the 
plateau, or peak value of G governs the delamination onset at the edge or, 
in this case, at the initial point of the taper [5], The peak values of G 
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on the right side in Figure 9 are plotted on the right side of Figure 10 as 
a function of the corresponding a/h. Similarly, the peak G values on the 
left side in Figure 9, occurring near point C, are plotted on the left side 
of Figure 10 as a function of the corresponding b/h 4 The numerical values 
are included in Table 4. The peak G values at a/h=0 (for b/h-4.5) and at 
b/h=0 (for a/h=5.75) are nearly equal and may be hypothesized as the 
critical value for the onset of a delamination at point C under fatigue 
loading. Delamination in the tapered laminate can be predicted by comparing 
this value to the threshold for delamination onset [6] . 

If the delamination initiates at point C under static tension loading, 
its growth will be governed by a mixed-mode criterion [5] because both the 
mode I and mode II components of G are present due to the tapered 
configuration. Figure 11 shows the percentages of mode I and mode II at 
delamination tip H in the thin laminate corresponding to a value of a-24h. 
Table 5 summarize results for several values of a/h. The mode I component is 
predominant for all values of b/h < 18. 

In contrast, as shown in Figure 12 and Table 6, the G components at 
delamination tip I in the tapered region CB initially consists of a large 
mode I component but at a/h=*18 the mode I component is only 50 percent of 
the total and continues to decrease with increasing a/h. Near the plydrops, 
Gj drops suddenly, but then increases. Overall, G^ decreases with increasing 

a/h, and G ^ increases with increasing, a/h. The value of G^ is 

approximately zero (i.e. G becomes 100 percent mode II) at a/h— 54. This 
distance depends on the initial delamination length, b/h, in the thin 
laminate . 

The mode I component of G Is predominant for a small delamination 
initiating at point C (Fig.l) and growing either along CD (Fig. 11) or along 
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CB (Fig. 12). The corresponding peak G^ values for various a/h and b/h ratios 

are shown in Table 7 and are plotted in Fig. 13. This figure is constructed 
in similar manner to Fig. 10. The peak values at a/h^O for growth into the 

thin region or at b/h~0 for growth into the tapered region may be compared 

to G r for the composite to predict delamination onset under static tension 
I C 

loading [ 5] . 


CONCEDING REMARKS 


A tapered composite laminate subjected to tension load was analyzed 
using the finite element method. The stacking sequence of the laminate was 
assumed to be { [ 0^/(±45) ] [ (±45) ^ ] /[0/(±45)/0] } g . The group of (±45)^ 

plies was dropped in three distinct steps, each 20 ply- thicknesses apart, 
thus forming a taper angle of 5.71 degrees. Neat resin pockets are assumed 
at the ends of ±45 degree plies that were terminated. The material of the 
laminate was assumed to be S2/SP250 glass/epoxy. 

A two-dimensional plane strain analysis was performed to determine 
stress distributions in the laminate without a delamination. The 
interlaminar normal stress and interlaminar shear stress distribution along 
the tapered interface, indicated by an arrow in the above stacking sequence, 
were calculated. Then delaminations were assumed to initiate at the point of 
intersection of the tapered interface and the thin region of the laminate. 
Delamination growth in the finite element model was simulated along the 
taper and into the thin region. The total strain-energy-release rate, G, and 


the mode I , 

and 

mode 

II components, G^ 

and Gjj, 

were computed at the 

delamination 

tip 

using 

the virtual crack 

closure 

technique (VCCT) . 
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Alternatively, G was obtained from a global energy balance. Based on the 
analysis performed here, the following conclusions were reached: 

1) Steep gradients of interlaminar normal and shear stress exist at the 
points of material and geometric discontinuities created by the internal 
plydrops. The largest value of interlaminar normal stress appears to occur 
at the intersection of the tapered interface and the thin region of the 
laminate. This is probably the site where a delamination would initiate. 

2) The strain-energy-release rate, G, was calculated for a delamination 
initiating at a point, located at the intersection of the taper and the thin 
laminate, and lying on the interface indicated by the arrow in the layup. 
The G values increase continually as the delamination grows into the thin 
laminate portion or along the taper. This indicates that a delamination 
initiating at the end of the taper will grow unstably along the taper and 
the thin laminate simultaneously. 

3) The strain-energy-release rate for a delamination growing a short 
distance into the thin laminate consists predominantly of mode I (opening) 
component . 

4) For a delamination growing along the tapered region, the strain- 
energy-release rate was initially all mode I but decreased with increasing 
delamination size until eventually it was all mode II. 

These results may help understand the delamination behavior in the 
tapered laminates and may be useful in predicting the onset and growth of 
the delamination under static and fatigue loading. 
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Appendix A 


Transformation of Stiffness Coefficients 
The stress -strain relations for each lamina were transformed from the 
material coordinate system 1,2,3 (Fig. 14) to the global system XYZ using 
the following procedure. The 3-D stress - strain relation for a ply in the 
material coordinate system is 

(a) 123 = [C] U) 123 (Al) 

where { cr } ^22 a 33 a \2 °23 a 13*’ 

(e) 123 “ l€ ll £ 22 e 33 € 12 € 23 £ 13 5 ’ 

and [Cl, , is a matrix that can be determined from elastic constants. 

L J 6x6 

Following similar notations, the stress-strain relations for a lamina in the 
global system can be written as 

|cr, XYZ " [C] (£} XYZ (A2) 


The matrix [C] is obtained from matrix [C] by rotating the material system 

1,2,3 (Fig. 14) to the global coordinate system XYZ through two rotations; 
a rotation ( 9 ) about the Z (or 3) axis, and then by a rotation ( <f > ) about 

r ' 

the Y (or Y ) axis. The transformed stiffness coefficient matrix, [C] , is 
obtained from the material stiffness coefficient matrix, [C] as 


[C] 


6x6 


IV6x6 'V6x6 t C '6x6 'VL 6 <VL 


(A3) 


where [T.] and [T,] are defined in terms of the appropriate angle as 
9 <p 
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2. 
cos 8 

sin^# 

0 

2*cos0*sin0 

0 

0 


2 

sin 6 

2 

cos $ 

0 

- 2*cos0*sin0 

0 

0 

[T ff ] - 

0 

0 

1.0 

0 

0 

0 


-cos0*sin0 

cos#*s in0 

0 

2 2 
cos 9 - sin 9 

0 

0 


0 

0 

0 

0 

cos0 -sin0 

and 

0 

0 

0 

0 

sin# cos$ 


2, 
cos <p 

0 

sin^ 

0 

0 

2*cos<£*sin<£ 


0 

1.0 

0 

0 

0 

0 

tV = 

sin^0 

0 

2. 
cos <p 

0 

0 

-2*cos<£*sin.£ 


0 

0 

0 

cos^ 

sin<£ 

0 


0 

0 

0 

-sin<£ 

cos<£ 

0 


-cos^*sin^ 

0 

cos^*sin<£ 

0 

0 

2 2 
cos <f> -sin <f> 


The superscript T to a matrix in equation (A3) denotes the transpose of the 
matrix. Furthermore, the plane strain conditions require that €yy“ e xy” € yz~ 
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0. Incorporating these conditions in (A2) yields the stress-strain relations 
for plane strain as 



{<7, xz ‘ 

^xz 

(£) xz 


(A4) 

where U) xz “ ^xx a zz a xz }; 

(£, xz 

■ {£ xx 

£ zz e XZ ) ' 


t 

and obtained from the 

global 

t 

[C] matrix as, 









c n 

C 13 

C 16 



t c lxz - 

C 31 

C 33 

C 36 


(A5) 


C 61 

C 63 

C 66 
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TABLE 1: 


MATERIAL PROPERTIES USED IN THE ANALYSES 


S2/SP250 GLASS/EPOXY : RESIN : 


E 11 


- 7.30 MSI 

E 

- 0.59 MSI 

E 22 ’ 

E 33 

- 2.10 MSI 

G 

- 0.224 MSI 

G 12 ’ 

G 13 ’ 

G 23 - 0.88 MSI 

V 

- 0.33 

l/ 12’ 

"l3’ 

!/ 23 - 0.275 




0.275 


TABLE 2: TOTAL STRAIN- ENERGY-RELEASE RATE AT DELAMINATION TIP 'H' 


IN THE THIN LAMINATE 


Gh/N 2 *10 12 


b/h 


a/h=0 a/h=6 a/h=12 a/h=>=20 a/h— 24 



.43 

10, 

.57 

110 

.00 

171 

.00 

230 

.71 

287 

.86 


.90 

17. 

.50 

113 

.58 

178 

.67 

242 

.83 

302 

.42 

1 

.50 

24. 

.17 

115, 

.75 

185 

.08 

254, 

.08 

316 

.17 

2. 

.10 

28. 

,42 

116, 

.08 

189 

.00 

262, 

.42 

326 

.67 

3. 

.00 

31. 

,79 

116. 

.21 

194 

.63 

275. 

.33 

343 

.54 

4. 

.20 

32, 

,79 

114. 

.33 

199, 

.00 

288. 

.83 

362 

.08 

5, 

.40 

31. 

,67 

111. 

,63 

201, 

.42 

299. 

.92 

378. 

.25 

6, 

.90 

29. 

,25 

107. 

.78 

203. 

.00 

311. 

81 

396. 

.61 

8. 

.70 

26. 

.08 

102. 

,81 

202. 

,83 

322. 

64 

414. 

.75 

10. 

.80 

22. 

.85 

96. 

,98 

200. 

.69 

331. 

40 

431. 

.29 

15. 

.00 

18. 

.15 

86. 

,74 

194. 

.06 

341. 

.78 

455. 

.43 

21. 

.00 

13. 

.59 

73. 

.30 

178, 

,81 

339. 

28 

465. 

.62 

30. 

,00 

12 . 

20 

64. 

98 

166. 

.61 

336. 

18 

474. 

.39 
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TABLE 3: 


TOTAL STRAIN- ENERGY-RELEASE RATE AT DELAMINATION TIP 


ALONG THE TAPER 


Gh/N 2 *10 12 


a/h 


b/h=0 b/h=6 b/h-12 b/h=24 



.75 

19 

.40 

93 

. 10 

173 

.00 

260 

.90 

2 

.75 

42 

.23 

130 

.46 

207 

.04 

333 

.14 

5 

.75 

47 

.82 

131 

.12 

224 

.82 

400 

.60 

7 

.50 

45 

.60 

126 

.15 

224 

.90 

424 

.00 

9 

.00 

43 

.50 

121 

.25 

222 

.48 

437 

.50 

11 

.00 

40 

.93 

113 

.55 

214 

.75 

445 

.25 

13 

.00 

38 

.58 

104, 

.63 

202 

.65 

440 

.25 

15 

.00 

37, 

.05 

95, 

, 13 

186 

.20 

421 

.25 

17 

.00 

42. 

.55 

94. 

55 

177, 

.88 

401. 

.50 

18. 

.50 

70, 

,30 

128. 

40 

214. 

.25 

439. 

,50 

19. 

.13 

115. 

40 

183. 

40 

276. 

,20 

510. 

00 

19. 

.38 

113. 

20 

181. 

40 

273. 

40 

502. 

00 

19. 

63 

132. 

60 

205. 

40 

300. 

20 

530. 

00 

19. 

88 

160. 

40 

238. 

60 

336. 

80 

570. 

00 

20. 

13 

275. 

20 

413. 

40 

590. 

00 

1004. 

00 

20. 

38 

298. 

20 

418. 

20 

577. 

00 

970. 

00 

20. 

63 

288. 

00 

397. 

20 

544. 

40 

916. 

00 
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TABLE 3 (contd.) 


2 12 
Gh/rr *10 


in. 

1ST 


2 


a/h 

b/h-0 

b/h=6 

b/h-12 

b/h-24 

21.00 

280.10 

383.50 

523.70 

882.00 

22.13 


336.20 

471.11 

822.86 

24.50 


256.55 

383.63 

729.75 

27.00 


207.35 

323.40 

657.03 

29.00 


184.30 

291.20 

611.98 

31.00 


167.33 

264.48 

569.05 

34.00 


148.30 

228.23 

497.18 

37.75 



241.60 

469.31 

39.63 



340.20 

578.00 

39.88 



380.20 

629.40 

40.25 



636.40 

1082.60 

40.75 



624.50 

1011.20 

42.00 



581.33 

935.78 

44.50 



483.47 

808.02 

48.00 




671.81 

51.00 




573.73 

53.00 




512.15 
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TABLE 4 : 

PEAK 

VALUES OF 

TOTAL STRAIN- 

ENERGY-RELEASE 

RATE 

DELAMINATION 

TIP 

a/h 

b/h 

PEAK G*h*10 12 

m 2 



N x 

lb. 

H 


0.0 

4.5 

32.5 


H 


6.0 

3.0 

116.5 


H 


12.0 

7.0 

204.0 


H 


20.0 

14.5 

340.0 


i 


5.75 

0.0 

47.5 


i 


3.00 

6.0 

130.0 


i 


6.00 

12.0 

220.0 


i 


11.00 

24.0 

450.0 



26 



TABLE 5a: PERCENTAGE OF Gj/G AT DELAMINATION TIP 'H' 

IN THE THIN LAMINATE 


b/h 


G. 

G 

I 

-- x 100 



a/h=0 

a/h=6 

a/h=12 

a/h-20 

a/h-24 

1.20 

99.05 

98.33 

97.69 

96.92 

95.74 

1.80 

97.69 

99.66 

99.19 

98.64 

97.73 

3.60 

94.57 

99.88 

100.00 

100.00 

99.88 

4.80 

92.16 

98.72 

99.50 

99.86 

100.00 

7.80 

86.53 

95.25 

97.07 

98.00 

98.91 

18.00 

73.01 

87.43 

91.40 

93.17 

94.89 
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TABLE 5b: PERCENTAGE OF G n /G AT DELAMINATION TIP ' H 1 


IN THE THIN LAMINATE 


G 


II 


x 100 


G 


b/h 


a/h=0 a/h=6 a/h=12 a/h=20 a/h=24 


1 . 

.20 

0 . 

.95 

1 . 

.67 

2 

.31 

3 

.08 

4 

.26 

1 . 

00 

o 

2. 

. 31 

0 . 

.34 

0 . 

.81 

1 

.36 

2. 

.27 

3. 

.60 

5. 

.43 

0 . 

, 12 

0 . 

.00 

0 

.00 

0. 

.12 

4, 

oo 

o 

7. 

, 84 

1 . 

,28 

0 . 

.50 

0 

.14 

0 . 

.00 

7. 

.80 

13. 

.47 

4. 

,75 

2. 

.93 

2, 

.00 

1 . 

.09 

18. 

o 

o 

26. 

.99 

12. 

57 

8. 

,60 

6. 

.83 

5. 

11 


TABLE 6a: PERCENTAGE OF G /G AT DELAMINATION TIP ' I ' 


ALONG THE TAPER 


a/h 


G 


I 


x 100 


G 


b/h =0 b/h -6 b/h =*12 b/h =24 



.50 

88 

.46 

99 

.25 

99 

.38 

99 

.20 

1 

.00 

86 

.15 

97 

.13 

97 

.80 

97 

.80 

4 

.50 

71 

.85 

83 

.68 

87 

.53 

89 

.70 

7 

.00 

56 

.77 

73 

.17 

79 

.66 

83 

.80 

10 

.00 

39 

.98 

62 

.67 

72 

.46 

79 

.10 

12 

.00 

29 

.09 

55 

.29 

67 

.48 

76 

.10 

14 

.00 

18 

.82 

47 

.35 

62 . 

.14 

73 . 

.00 

16 

.00 

10 , 

.09 

37 . 

.74 

55 , 

,11 

68 . 

.90 

19 . 

,25 

11 . 

.61 

26 . 

,88 

40 . 

54 

56 . 

00 

19 . 

50 

13 . 

27 

27 . 

33 

40 . 

01 

54 . 

80 

19 . 

75 

14 . 

27 

26 . 

50 

37 . 

67 

51 . 

30 

20 . 

00 

13 . 

95 

25 . 

09 

35 . 

10 

47 . 

50 

20 . 

25 

6 . 

61 

18 . 

69 

30 . 

59 

45 . 

80 

20 . 

50 

2 . 

00 

11 . 

00 

22 . 

26 

38 . 

60 

20 . 

75 

1 . 

01 

8 . 

70 

19 . 

53 

36 . 

20 

21 . 

25 


54 

7 . 

36 

17 . 

91 

34 . 

70 
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TABLE 6b: PERCENTAGE OF G n /G AT DELAMINATION TIP 'I' 


ALONG THE TAPER 


a/h 


G II 

G 

x 100 


b/h-0 

b/h=6 

b/h=12 

b/h-24 

.50 

11.54 

.75 

.62 

.80 

1.00 

13.85 

2.87 

2.20 

2.20 

4.50 

28.15 

16.32 

12.47 

10.30 

7.00 

43.23 

26.83 

20.34 

16.20 

10.00 

60.02 

37.33 

27.54 

20.90 

12.00 

70.91 

44.71 

32.52 

23.90 

14.00 

81.18 

52.65 

37.86 

27.00 

16.00 

89.91 

62.26 

44.89 

31.10 

19.25 

88.39 

73.12 

59.46 

44.00 

19.50 

86.73 

72.67 

59.99 

45.20 

19.75 

85.73 

73.50 

62.33 

48.70 

20.00 

86.05 

74.91 

64.90 

52.50 

20.25 

93.39 

81.31 

69.41 

54.20 

20.50 

98.00 

89.00 

77.74 

61.40 

20.75 

98.99 

91.30 

80.47 

63.80 

21.25 

99.46 

92.64 

82.09 

65.30 
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TABLE 6b (contd.) 


G 


II 


- x 100 


G 


ct/ U - - 

b/h=0 b/h=6 

b/h-12 

b/h-24 

23.00 

90.41 

78.49 

60.70 

26.00 

89.67 

76.16 

57.00 

28.00 

91.33 

77.45 

57.20 

30.00 

93 . 64 

79.83 

58.30 

32.00 

96.15 

82.98 

60.20 

36.00 

99.91 

91.33 

67.40 

39.50 


90.39 

73.20 

39.75 


89.42 

73.90 

40.00 


86.89 

72.30 

40.50 


93.72 

78.80 

41.00 


96.64 

83.40 

43.00 


97.25 

83.70 

46.00 


98.24 

84.00 

50.00 



89.80 

52.00 



94.30 

54.00 



98.90 
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TABLE 7: 


PEAK VALUES OF mode I STRAIN-ENERGY-RELEASE RATE 


DELAMINATION TIP 

a/h 

b/h 

Gjh 

* io 12 in . 2 /l 

N 2 

X 

H 

0.0 

3.60 

30.12 

H 

6.0 

1.80 

115.67 

H 

12.0 

4.80 

198.68 

H 

20.0 

7.80 

309.81 

I 

4.5 

0.0 

35.02 

I 

4.5 

6.0 

111.75 

I 

4.5 

12.0 

194.11 

I 

7.0 

24.0 

345.76 
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Fig. 2(b): Mesh detail showing delamination tip at point 
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stress distribution along 
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Fig. 6: Details of geometric and material discontinuities along 
interface ABCD 
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Tapered plane 'BC\ delam. tip "I" a Plane 'CD', delam. tip "H 

Point C 

Fig. 9: Normalized total strain-energy-release rate at delamination 
tip I along interface CB and tip H along interface CD 



450 x 10 
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Fig. 10: Peak values of total strain-energy-release rate along either 



Fig. 11: Gj/G and G^/G at delaraination tip H along interface CD 
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Fig. 12: G t /G and G TT /G at delamination tip I along interface BC 




Fig. 13: Peak values of G, on either side of C 




Fig. 14: Material and global coordinate systems and rotations 
used in the transformations 


NASA 

Naiooa» Ae^orauics aro 
Soace Ao^^siraicri 


Report Documentation Page 


1. Report No. 

NASA m- 102592 
AVSCOM TM-90-B-003 


2. Government Accession No. 


4. Title and Subtitle 


Strain Energy Release Rate Analysis of Delamination 
in a Tapered Laminate Subjected to Tension Load 


7. Authorls) 


S. A. Salpekar*; I. S. Raju**; and T. K. O'Brien*** 


9. Performing Organization Name and Address 

NASA Langley Research Center, Hampton, VA 23665-5225 
U.S. Army Aviation Research Technology Activity (AVSCOM) 
Aerostructures Directorate 
Hampton, VA 23665-5225 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 

Washington, DC 20546 

U.S. Army Aviation Systems Command 

St. Louis, MO 63166 


3. Recipient's Catalog No. 


5. Report Date 

March 1990 


6. Performing Organization Code 


8. Performing Organization Report No. 


10. Work Unit No. 

505-63-01-05 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 


(r c t~- 45*. 


*S. A. Salpekar, Analytical Services and Materials, Inc., Hampton, VA 23666 

**I. S. Raju, North Carolina A & T State University, Greensboro, NC 27411 

***T. K. O’Brien, U.S. Army Aviation Research and Technology Activity (AVSCOM), Langley 

Research Center, Ha mpton, VA 23665 

16. Abstract 

A tapered composite laminate subjected to tension load was analyzed using the finite-element method. The 

{ [07/(±45)]/ 1 [(±45)3]/[0/(±45)/0]} s glass/epoxy laminate has a (+45)3 g rou P of plies dropped in three 

distinct steps, each 20 ply-thicknesses apart, thus forming a taper angle of 5.71 degrees. Steep gradients 
of interlaminar normal and shear stress on a potential delamination interface suggest the existence of stress 
singularities at the points of material and geometric discontinuities created by the internal plydrops. The 
delamination was assumed to initiate at the thin end of the taper on ffie -45/+45 interface indicated by the 
arrow in the laminate layup and the delamination growth was simulated in both directions, i.e., along the 
taper and into the thin region. The strain-energy-release rate for a delamination growing into the thin lami- 
nate consisted predominantly of mode I (opening) component. For a delamination growing along the 
tapered region, the strain-energy-release rate was initially all mode I, but the proportion of mode I 
decreased with increase in delamination size until eventually total G was all mode IT. The total G for 
both delamination tips increased with increase in delamination size, indicating that a delamination initiating 
at the end of the taper will grow unstably along the taper and into the thin laminate simultaneously. 


fafMlMi u ^r ,ed bv Au,ho,,s ” 

Plydrop 

Delamination 

18 Distribution Statement 

Unclassified - Unlimited 

Strain-energy-release rate 
Two dimensional analysis 

Subject Category - 24 


19. Security Classif. (of this reportl 20 Security Classif. (of this page) 21 No. of pages 22. Price 

Unclassified Unclassified 50 A0 3 


NASA FORM 1826 OCT 86 









